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Using a differential technique, we have measured the specific heats of polycrystalline 
Bai_ a: K a; Fe2As2 samples with x = 0, 0.1 and 0.3, between 2K and 380K and in magnetic fields 
to 13 Tesla. From this data we have determined the electronic specific heat coefficient 7(= C e i/T) 
over the entire range for the three samples. The most heavily doped sample (x — 0.3) exhibits a 
large superconducting anomaly A"/(T C ) ~ 48 mj/mol K 2 at T c — 35K, and we determine the energy 
gap, condensation energy, superfluid density and coherence length. In the normal state for the x = 
0.3 sample, 7 ~ 47 mj/mol K 2 is constant from T c to 380K. In the parent compound (x = 0) there 
is a large almost first order anomaly at the spin density wave (SDW) transition at T = 136K. This 
anomaly is smaller and broader for x = 0.1. At low T, 7 is strongly reduced by the SDW gap for 
both x — and 0.1, but above T a , 7 for all three samples are similar. 



PACS numbers: 74.25.Bt, 74.70.-b 

The electronic specific heat contains a wealth of quan- 
titative information on the electronic spectrum over an 
energy region zblOOmeV about the Fermi level, cru- 
cial to the understanding of high temperature supercon- 
ductivity. Measurements of the electronic specific heat 
played an important role in revealing key properties of 
the copper-oxide based 'cuprate' high-temperature su- 
perconductors (HTSCs). Some examples include the 
normal-state "pseudogap" [IHlj , the bulk sample inho- 
mogeneity length scale [5], and more recently the degree 
to which the superconducting transition temperature is 
suppressed due to superconducting fluctuations |B]. In 
this work we extend such measurements to the recently 
discovered iron-arsenide based 'pnictide' HTSCs. Here 
we present results obtained from polycrystalline samples 
of Bai_ a; K :c Fe2As2 (x = 0, 0.1 and 0.3) using a high- 
resolution differential technique [7 . 

With this technique we directly measure the difference 
in the specific heats of a doped sample and the undoped 
reference sample (BaFe2As2). This eliminates most of 
the large phonon term from the raw data and yields a 
curve which is dominated by the difference in electronic 
terms between the sample and reference. After making a 
small correction for any residual phonon term, this differ- 
ence in electronic specific heats can be determined with 
a resolution of ~ 0.1 mj/mol K 2 at temperatures from 
1.8K to 380K and magnetic fields from to 13T. During 
a measurement run the total specific heats of the sample 
and reference are also measured. 

The polycrystalline samples of Bai_ 3; K 2 ;Fe2As2 were 
prepared by a solid state reaction method similar to that 
reported by Chen et al 8\. First, Fe 2 As, BaAs, and 
KAs were prepared from high purity As (99.999%), Fe 
(99.9%), Ba (99.9%) and K (99.95%) in evacuated quartz 
ampoules at 800, 650 and 500°C respectively. Next, 
the terminal compounds BaFe2As2 and KFe2As2 were 
synthesized at 950 and 700° C respectively, from stoi- 
chiometric amounts of BaAs or KAs and Fe2As in alu- 



mina crucibles sealed in evacuated quartz ampoules. Fi- 
nally, samples of Bai_ a; K. E Fe2As2 with x = 0.1 and 0.3 
were prepared from appropriate amounts of single-phase 
BaFe2As2 and KFe2As2- The components were mixed, 
pressed into pellets, placed into alumina crucibles and 
sealed in evacuated quartz tubes. The samples were an- 
nealed for 50 h at 700° C with one intermittent grind- 
ing, and were characterized by room temperature pow- 
der X-ray diffraction using Cu K Q radiation. The diffrac- 
tion patterns were indexed on the basis of the tetragonal 
ThCr2Si2 type structure (space group I4/mmm). Lattice 
parameters calculated by a least-squares method agree 
well with those reported by Chen et aZ[5]. Traces of FeAs 
as an impurity were detected for compositions with x — 
0.1 and 0.3. The samples for heat capacity measurement 
weighed ~ 0.8g. 

The total specific heats of the three samples are plot- 
ted as 7 to * = C tot /T in Fig. [TJ In the x = sample 
there is a sharp and almost first order anomaly at the 
magneto-structural transition at To = 136K, in agree- 
ment with published single crystal data|9j. For the x = 
0.1 sample the corresponding anomaly at T = 135K is 
broader and considerably reduced in magnitude, in agree- 
ment with the data of Rotter et al. [TU] who have tracked 
To in the specific heat down to 105K in a sample with x 
= 0.2. The magnetic field dependence of this anomaly is 
extremely weak for both samples. We see no evidence for 
a structural transition in the x = 0.3 sample (the very 
weak anomaly at 67K in the differential data (Fig. [2]ja)) 
is probably due to an FeAs impurity phaseJT]). The su- 
perconducting transition in this sample is clearly visible 
in 7*°* at 35K, and the field dependence is shown in the 
inset to Fig. [T] 

In the temperature range 2 - 8K, j tot = 7 (0) + /3(0)T 2 
at zero field with 7(0) = 5.1, 9.3 and 1.8 mj/mol K 2 and 
/3(0) = 0.45, 0.37 and 0.57 mj/mol K 4 for the x = 0, 0.1 
and 0.3 samples respectively. On the assumption that 
/3(0) is entirely due to phonons we obtain Debye temper- 
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FIG. 1: (Color online) 7'°* = C tot /T of Bai-^FeaAsa for 
x = 0, 0.1 and 0.3. The dotted line shows the high temper- 
ature limiting value of the lattice contribution, 3nR/T. The 
inset shows the suppression of the superconducting anomaly 
in magnetic fields from 1 to 13T in IT increments. 



atures 0^(0) = 289, 298 and 257K for the three samples. 
We cannot, however, ignore the possibility that for the 
superconducting 0.3 sample /3(0) is enhanced by an elec- 
tronic component oc T 2 , leading to too low a value of 
#d(0). This would explain the reduced value of #d(0) 
for this sample which contradicts the trend seen in A7*°* 
(Fig.[2ja)) that acoustic phonon frequencies, which dom- 
inate the phonon specific heat in this temperature region, 
increase with doping. If the and 0.1 samples have an 
initial electronic T 2 term in "f(T) we would expect this 
to have a magnitude ~ 7„T 2 /T 2 ~ 0.003 mJ/mol K 2 
if it is controlled by the same energy scale as To. This 
contribution to /3 would be negligible compared with the 
measured /3(0). 

We next discuss the separation of electronic and 
phonon contributions from the total specific heat over 
a wider temperature range. The total specific heat co- 



efficient is given by 7 



7 + 7 p,i + 7°™ where the 
electronic term 7 = C e ;/T, the harmonic phonon term 
j' ph = C v /T and the anharmonic phonon (dilation) term 
7 Q « = (C p - C V )/T]T%. C p and C v are the lattice 
heat capacities at constant pressure and volume respec- 
tively. The anharmonic term is given by 7°™ = VBB 2 ^T7\ 
where V is the molar volume, B is the bulk modu- 
lus and f3 is the volume expansion coefficient. We as- 
sume that 7°" is doping independent and use V = 61 
cm 3 /mol[in], B ~ 0.80 xlO 8 mJ/cm 3 13], and /3(300K) 
~ 50xl0 -6 /K[T1] to obtain a room temperature value 
7 an (300K) - 12 mJ/mol K 2 for each sample. To a 
very good approximation (3(T) cx C V (T) 12J and thus 
7 Q "(T) = [C V (T)/C V (300K)] 2 ■ 7 an (300iT). A plot of 
7 a "(T) is shown in the inset to Fig. [2|a). 

After correcting for j an (T), the electronic term 7 can 
be determined at high temperatures where C v is close 



to its saturation value 3ni? (where n is the number of 
atoms/formula unit) shown by the dotted line in Fig. [I] 
To extend the useful range of the high temperature region 
we exploit the fact that the Debye temperature 0d{T) 
deduced from C v is generally only weakly T-dependent 
for T > 0.5#£i(oo) (or ~ 180K for the present materials). 
With this constraint we obtain limiting high temperature 
Debye temperatures #d(oo) = 368, 360 and 356K and 
high-T values 7(300^) = 55, 47 and 47 mJ/mol K 2 for 
the 0, 0.1 and 0.3 samples. 

Differential measurements (see Fig. ^ a)) between each 
sample and the x = reference give Aj tot — A7 + Aj?* 1 
(assuming A~/ an = 0). It is clear from Figs. [TJ and 
[5Ja) that the phonon terms for the three samples are 
very similar. The broad negative peak at 35K seen in 
A7* ot (0.1,0) = j tot (x = 0.1)-7*°*(x = 0), and which can 
also be inferred in the data for A7* ot (0.3,0) in the same 
temperature region, is compatible with fractional in- 
creases of 0.018 and 0.030 in the acoustic phonon frequen- 
cies for the 0.1 and 0.3 samples relative to the x = ref- 
erence. Since we expect the difference of electronic terms 
A7(0.1,0) for the two non-superconducting samples to 
be only weakly T-dependent (at least below 100K), we 
obtain A 7 p /i (0.1,0) = A7*°*(0.1,0) - A 7 (0.1,0) using 
the low temperature value A7(0.1,0) = —4.2 mJ/mol 
K 2 . To ensure that the resulting A7 ph (0.1,0) has a 
T-dependence compatible with that of a phonon spec- 
trum it was modelled with a histogram for the difference 
phonon spectrum, as discussed previously |15j. Making 
this phonon correction removes the broad negative peak 
in A7* ot at 35K for this sample (Fig. [2](b)) . As shown in 
the same figure, subtracting a phonon term A7 pft- (0.3,0) 
= 1.68A7 p/l (0.1,0) also removes the negative peak in 
A-f tot (0.3,0). Following this correction, the electronic 
A7(0.3,0) has an additional negative T-dependence given 

by 20 x 10~ 6 T 3 mJ/mol K 2 in the range 40 - 110K. 

We conclude below that 7«(T) for the 0.3 sample is T- 
independent, and attribute this negative T-dependcncc 
to a corresponding positive term in the electronic spe- 
cific heat coefficient 7(2; = 0) of the undoped sample in 
the SDW phase. 

After correcting for A7 p,i (0.3,0) we find that just above 
T c , j(x = 0.3) - 47 mJ/mol K 2 . This is very close to the 
value S/T(T C ), where S = J Q T "f(T)dT is the electronic 
entropy. By definition, S/T(T') is the average value of 7 
below T' . From conservation of entropy S/T(T C ) is also 
the average value of the underlying normal state 7n(T) 
below T c . Therefore the result above is consistent with 
a T-indepcndent 7 ?l (T) for T < T c . In addition, this 
value for 7«(T) below T c is close to its high temperature 
value estimated above, and it is reasonable to assume 
that for this sample 7n(T) is T-indepcndent over the en- 
tire temperature range. Subtracting this constant term 
from j tot (x = 0.3) gives j ph (x = 0.3), which we then fit 
with a nine-term histogram for the phonon spectrum|15j. 
Finally, subtracting A 7 p /i (0.3,0) and A 7 p' 1 (0.3,0.1) from 
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FIG. 2: (Color online) (a) The directly measured raw differ- 
ence in specific heats between the doped samples and undoped 
reference in and 13T fields. The inset shows the anhar- 
monic phonon term 7 an . (b) The electronic specific heat of 
Bai_ a; K ;c Fe2As2 for x = 0, 0.1 and 0.3. 

j ph (x = 0.3) gives j ph (x = 0) and ^(x = 0.1). 

The electronic specific heats of the three samples are 
shown in Figure [2jb). Above 150K, 7 is at most only 
weakly T-dependent and almost independent of doping. 
At low temperatures 7 for the x = and 0.1 samples 
is heavily reduced by factors of ~ 11 and 5 respectively 
due to the opening of a spin density wave (SDW) gap be- 
low T ~ 136K, and increases as - 20 x 10" 6 T 3 mJ/mol 
K 2 in the range 40 — 110K. Band splitting and signs of 
partial gapping of the Fermi surface have been observed 
in the SDW state of BaFe2As2 by angle-resolved photoc- 
mission spectroscopy|16] . 7 for the x = 0.3 sample is 
dominated by the large superconducting anomaly with 
A7(T C ) ps 48 mJ/mol K 2 , a normal-state 7at(0) of about 
47 mJ/mol K 2 and 7(0) = 1.8 mJ/mol K 2 . Measure- 
ments on a more overdoped crystal with x = 0.4 by Mu 
et aZ.[T7] show an even larger Aj(T c ) 100 mJ/mol K 2 
and 7at(0) f=s 63 mJ/mol K 2 , implying a further growth 
in the density of states with doping. For our x =0.3 
sample the superconducting condensation energy U (0) = 
10.6 J/mol, as determined from J Q " (S n — S s ) dT, where 
S n and S s are the normal and superconducting state en- 
tropies respectively. In Fig. [2jb), the curvature in 7 for 
x =0.3 between 10 and 20K is reduced compared to the 
raw Aj tot data (Fig. [2]ja)) by the phonon correction de- 
scribed above. Evidence that the remaining curvature 
in 7 may reflect a genuine anomalous T-dependence of 
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FIG. 3: (Color online) Temperature dependence of the super- 
conducting gap extracted by comparing the experimentally 
determined S/T with that calculated from models of the den- 
sity of states (shown in the inset) given by Eqns. [TJ3] 

the electronic term comes from the observation of similar 
anomalous curvature between 15 and 20K in the field de- 
pendence Aj(H) — -f(H, T) — 7(0, T) shown in the inset 
to Fig. [4]j a). Since the phonon term is -ff-independent, 
this anomalous curvature in A'y(H) can only be of elec- 
tronic origin, and may signal the presence of a second 
energy gap as inferred by Mu et al[17\. 

To determine the magnitude and temperature depen- 
dence of the superconducting gap A for x = 0.3 (shown 
in Fig. [3]) we first subtract from 7 the small residual 7(0) 
= 1.8 mJ/mol K 2 and then match the corrected S/T 
with that calculated assuming the following models for 
the density of states (see Fig. [3] inset) : 

s-wave 

N(E) = E/^E 2 - A 2 (1) 

d-wave 

N(E) = ^rJ2 E/VE 2 - A 2 cos 2 29 (2) 
6 e 

and s + id-wave 

N(E) = j^^2 E /^ E2 - 0-98A 2 cos 2 26 - 0.02A 2 (3) 
6 9 

where the summations over 9 are from to 45° and Ng 
is the number of 9 values. Although such modelling is 
too crude to pin down the exact nature of the gap, the 
rapid increase in A below 15K in the case of the as- 
sumed d-wave gap function, shows that nodes on the 
Fermi surface are incompatible with the T-dependcncc 
of 7 observed at low temperature, even in the presence 
of weak pair breaking. The s + id- wave gap function pro- 
duces a more reasonable A(T). This function completely 
gaps the Fermi surface but is strongly anisotropic. Such 
anisotropy could arise from the presence of multiple gaps 
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FIG. 4: (Color online) (a) Magnetic field dependent change 
in electronic specific heat j(H) — 7(0) vs temperature, in IT 
increments, from the data in the inset of Fig. [I] (Inset) En- 
largement of the to 30K region, (b) Field dependence of 
the electronic entropy, S(H) — S'(O), calculated from (a) using 
Eqn.[5] (c) Field dependence of the free energy, F(H) — F(0), 
calculated from (b) using Eqn.|6] (Inset) A fit to F(H)-F(0) 
at 2.5K using Eqn. [|] 



on different sheets of the Fermi surface, or from an in- 
trinsically anisotropic gap(s). 

We turn now to the information contained in the mag- 
netic field dependence of 7. The inset to Fig. [T] shows 
the suppression of the superconducting anomaly in mag- 
netic fields from 1 to 13T in IT increments. Because the 
phonon specific heat is independent of magnetic field, the 
change in electronic specific heat with field is obtained by 
subtracting the zero field data from the data measured 
in a field (see Fig. |4]^a)). 



A 1 {H,T)= 1 {H,T)- 1 {Q,T) 



(4) 
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FIG. 5: Pollycrystalline averaged mixed state parameters of 
Bao.3Ko.3Fe2As2 derived from the field dependence of the 
free energy AF(H,T): (a) Critical fields H c , H c i, H C 2, and 
Ginsburg-Landau parameter k; (b) London penetration depth 
A, superfluid density oc 1/A 2 and superconducting coherence 
length £. 
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tronic entropy with field (see Fig.Qb)). 

AS(H,T) = [ A-f(H, T)dT 
Jo 

The change in free energy is obtained by integrating the 
entropy over temperature (see Fig.Qc)). 

AF(H,T) = - [ AS(H,T)dT + AF{H,0) (6) 
^0 

AF(H, 0) is determined from the condition that the su- 
perconducting contribution to AF(H,T) tends to zero 
for T T c . At each temperature, the field dependence 
of the free energy is fitted to a theoretical expression de- 
rived from the model of Hao and Clem[TS] for an s-wave 
superconductor. 



AF SjT . et ,(_ff, T) — 



a4>a 
32tt 2 A 2 



ef3H, 
H 



<2 



a x H\sx{a 2 H) 



(7) 
(8) 



Integrating over temperature gives the change in elec- 



The coefficients a and /3 are weakly field dependent and 
are given by a ~ 0.77 and j3 ~ 1.44 in the range 0.02 < 
H/H C 2 < 0.3. A fit to the free energy at 2.5K is shown 
in the inset to Fig. Qc). 

The penetration depth, A, and the upper critical field, 
H c2 , are determined directly from the fit parameters a\ 
and 02 . Then using the following Ginsburg-Landau rela- 
tions we extract the: critical field — (<^o/47rA 2 )_ff C 2, 
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Ginsburg-Landau parameter k = H c ij 'H c y2, lower crit- 
ical field H c i = (0o mK )/(47rA 2 ), and superconducting 



coherence length £ 



b /2TTH c: 



a/2 



The temperature 



dependencies of these quantities is plotted in Fig. [5] The 
values shown are polycrystalline averages. The large k « 
f 30 indicates the strong type II nature of this material. 
Our H C 2 values are very similar to those obtained from 
radio frequency penetration depth measurements on a 
single crystal with x = 0.45 [19 , and the value we ob- 
tain for A(0) « 260nm is similar to those measured by 
infrared spectroscopy [2U] and tunnel diode resonator pTj 
techniques. 

In summary, above 150K the electronic specific heats 
for the three samples are large and almost identical. For 
the x — sample there is a large and almost first order 
anomaly at the SDW transition at To = 136K, whilst the 
corresponding anomaly for the x = 0.1 sample at To = 
135K is smaller and broader and more closely resembles 
a second order phase transition. At low temperatures, 7 
for these two samples is strongly reduced by the SDW 
gap by factors of 11 and 5 for x = and 0.1 respectively. 
In the x = 0.3 sample the large superconducting anomaly 
and associated value of S/T(T C ) shows that the under- 
lying normal state 7 at low temperature is close to its 
high temperature value. This suggests that as x reduces, 
the growth of the SDW gap and consequent reduction in 
the density of states is responsible for the disappearance 
of superconductivity. The temperature dependence of 7 
(x = 0.3) supports a nodeless superconducting gap func- 
tion that is strongly anisotropic about the Fermi surface, 
possibly due to the presence of multiple gaps. 
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